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Abstract 

In this article we report on our results about quark produetion and ehemical equilibration of quark- 
gluon plasma. Our initial eondition eorresponds to a classie Yang-Mills spectrum, in which only 
gluon degrees of freedom are eonsidered; the initial eondition is then evolved to a quark-gluon 
plasma by means of relativistie transport theory with inelastic processes which permit the eon- 
version of gluons to qq pairs. We then compare our results to the ones obtained with a standard 
Glauber model initialization. We find that regardless of the initial eondition the final stage of the 
system eontains an abundant pereentage of qq pairs; moreover spanning the possible eoupling 
from weak to strong we find that unless the eoupling is unrealistically small, both production rate 
and final percentage of fermions is quite large. 
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1. Introduction 


In the last decade it has been reaehed a general consensus that Ultra-relativistie heavy-ion col¬ 
lisions (uRHICs) at the Relativistie Heavy-Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) ereate a hot and dense strongly interaeting quark and gluon plasma (QGP) [[ll,0,[E0]- A 
main diseovery has been that the QGP has a very small shear viscosity to density entropy, q/s, 
whieh is more than one order of magnitude smaller than the one of water [11,0], and close to the 
lower bound of l/4;7r eonjeetured for systems at infinite strong eoupling |13]. Aeeording to the stan¬ 
dard picture of ultrarelativistic heavy ion collisions, before the eollision the two colliding nuclei 
ean be represented as two thin sheets of eolor-glass condensate liiH whieh produee, immedi¬ 
ately after the eollision, a configuration of strong longitudinal chromoelectric and chromomagnetic 
fields named the glasma, see lilL O, [ij] for reviews. 

An interesting problem of uRHICs is the dynamical evolution of the high energy system made 
mainly of gluons, obtained from the deeay of the glasma flux tubes, to a locally thermalized 
and eventually chemically equilibrated quark-gluon plasma. This problem has been diseussed 
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previously in rtl4l] where quark-antiquark produetion rate is eomputed by means of the solution of 
the Dirae equation in the baekground of the strong initial glasma field, and in 111511 by simulations 
based on relativistie transport theory (RTT) whieh is a fruitful theoretieal tool to study the evolution 
of QGP produced in heavy ion collisions [16,17,18,1^^21,22,15, 23,1^. In this article 
we follow closely Ref. lilSn and in some sense the present study is a continuation of [jlSll : in 
fact in [jlSil several aspects have not been investigated, in particular the role of the initial non¬ 
equilibrium distribution on qq production times and on the chemical equilibration of the QGP, as 
well as a study of the coupling dependence of the aforementioned quantities. These are the aspects 
we study in the present article. 

We are interested to compute the qq production rate initializing simulations by a pure gluon 
plasma with the spectrum computed within the classical Yang-Mills (GYM) theory of the glasma. 
The evolution of the initial condition is then achieved by relativistic transport theory. The simula¬ 
tions with the GYM spectrum are started at tq = 0.2 fm/c, hence assuming in this time range the 
initial longitudinal gluon fields have decayed (the decay time is of the order oil/Qs with cor¬ 
responding to the saturation scale) and to populate the transverse momentum space (in the initial 
glasma the fields are purely longitudinal hence the pr-spectrum at r = 0"^ is zero). Within our 
approach we cannot discuss if and how quarks are produced before r = tq; however we find that 
even neglecting the possible quark formation for r < tq the QCD inelastic processes are efficient 
enough to obtain a final state which consists mainly of qq pairs rather than gluons. 

We also consider the problem of the chemical equilibration of the QGR For a system of quarks 
and gluons thermalized at the same temperature T the equilibrium value for R = {Nq + Nq)/Ng is 
given by 


R 


eq 
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( 1 ) 


for example if quarks and gluons are massless one finds easily R^q = 9/4 at chemical equilibrium. 
In the case of finite masses the ratio R depends on temperature even if quasiparticle masses are 
temperature independent. We find that changing the initial distribution from a Glauber model to 
a GYM spectrum has some effect on the chemical equilibration, the effect being more relevant in 
the case of weak coupling. 

In the present study we use the quasiparticles to identify quarks and gluons propagating in the 
QGR In particular we consider the case of massless quasiparticles, which corresponds to a fluid 
with the perfect massless gas equation of state; the case of massive quasiparticles with temperature 
independent masses, corresponding to a fluid with the equation of state of an ideal massive gas. 
Finally we also consider the case of energy density dependent quasiparticle masses, with masses 
fixed by requiring the equation of state of the fluid is the one of QCD as it is measured on the 
lattice 1 I 25 I] . We consider three quark flavors in our simulations. Finally, we make use of RTT with 
a Boltzmann kernel type with 2^2 processes, see Eq. @, hence neglecting the quantum nature of 
particles as well as bremsstrahlung processes. A comment on the latter follows: in the bottom-up 
scenario [l26n one would expect a significant increase of soft gluons due to gluon radiation; however 
it has been shown in ll27n that at RHIC and LHC conditions the direct and inverse processes nearly 
balance in the bulk, and the yield of gluons is not very affected by bremsstrahlung. Of course a 
quantitative study has to be done in order to give a firm statement about this important point, hence 
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we will address the problem of soft gluon production by bremsstrahlung in the near future. 

The plan of the article is as follows. In Section II we review relativistic transport theory which 
is the base of our simulation code. In Section III we summarize our results concerning formation 
time of quarks as well as chemical equilibration of the QGP, when we initialize simulations starting 
either from a Glauber model or a GYM spectrum; in this Section we consider the cases of massless 
as well as massive quasiparticles, in the latter case considering both constant and temperature 
dependent masses. Finally in Section IV we draw our conclusions. 


2. Relativistic Transport Theory 


In our study we simulate the evolution of the QGP fireball by means of a code based on the 
solution of the relativistic kinetic Boltzmann-Vlasov equation for the gluon and quark distribution 
function, 

[ifd, + M{x)d,M{xW] fix, p) = C[/] , (2) 

where = <9^, denotes the gradient with respect to momenta and Mix) corresponds to the quasi¬ 
particle mass: we have made the dependence of this mass on the space-time coordinates explicit, 
which comes from the fact that in quasiparticle models M is a function of energy density 11251] and 
the latter changes in space and time due to the expansion as well as to the interactions. It is useful 
to remind at this point that the use of a quasiparticle mass in the above equation implies that the 
QGP fluid has a given equation of state: if we take M = 0 it means we simulate a fluid with a per¬ 
fect gas equation of state s = 3P; on the other hand if M is fixed with a temperature dependence 
according to 11251] then it is equivalent to state that the QGP fluid evolves with the lattice QCD 
equation of state once the local kinetic equilibrium is reached. 

In Eq. Q the collision integral takes the form 


C[/] 
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(3) 


where A\ denotes the invariant amplitude for the process 12 ^ 1'2' and lEtdYi = d^pi/ilnf. 
Both elastic as well as inelastic processes contribute to the invariant amplitude. In this work we 
follow lll5l] and we consider only the two body processes. In particular all the elastic cross sections 
we take have the form 

dcXelastic ^ 

— T .— 7- 

dt it - mjff 

where t denotes the invariant Mandelstam variable and mo = 0.7 GeV is used as an infrared 
regulator and determines how isotropic the cross section is (the smaller the more forward 
peaked the cross section). The overall factor in the previous equation changes channel by channel 
iqq, qq, qg, qg and gg) as in UJ]. 

In this study we fix the value of and compute the relevant cross sections. In order to estimate 
an appropriate value of as we use the total cross section for gg scattering among massless gluons 
and require that the value of q/s obtained by this massless gluon-gluon scattering is approximately 
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the one used in hydrodynamies and/or transport simulations. To estimate sueh a eoupling we use 
the approximate kinetie theory relation ^3, 24] 
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whieh, for a reference temperature ofT = 0.3 GeV gives = 0.7 for 77/5 = IjAn as used in heavy 
ion collisions simulations 11161 l22L 12811 . A value of « 0.9 would be required for rjls = l/ 4 ; 7 r, 
which certainly would cause the gluons-to-quarks processes to be even faster but which does not 
affect considerably the final quark spectra as well as the final gluon number to quark number ratio, 
see next Section. We will also consider the case of a running strong coupling, 


aAT) = 


In 
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with A = 0.2 GeV; for example at T = 0.3 GeV one finds as ~ 0.4. 

For the inelastic processes between quarks and gluons gg <-4 qq we have A1 = Ms + Mt + Mu- 
For the massless case the cross sections for such processes are the textbook pQCD cross section 
for jet production in high-energy proton-proton collisions. With massive quarks the calculations 
are the Combridge cross sections |t29i|30l] used to evaluate heavy-quark production. In our case we 
have considered a finite mass for both gluons and quarks together with a dressed gluon propagator. 
It is worth to mention here that the cross section for the process at hand is dominated by the t- and 
M-channel and their interference, while the 5 -channel alone is negligible. The squared matrix 
element of the t-channel is given by 


|(t - [(mj - t)(mj -u)- 2ml(t + mj) - 4mJmJ - mj] , 

where mg, mg denote the quasiparticle masses for quarks and gluons respectively. 


( 7 ) 


3. Results 

3.1. The massless case 

In the present work we focus on two different initial conditions for our simulations. The first 
one is named thermal Glauber (Th-Glauber) in which we distribute particles in transverse coordi¬ 
nate space according to the usual Glauber model with the standard mixture O.SSNpan + O.lSVco//, 
while distributing particles in boost invatiant way along the longitudinal direction; initial momen¬ 
tum space distribution corresponds to a thermal distribution in transverse momentum space and a 
boost invariant distribution in the longitudinal direction with y = q, with y and q denoting momen¬ 
tum and space-time rapidities respectively. The other initial condition we use in our simulations is 
named GYM initialization, in which transverse plane particle distribution corresponds to a Glauber 
model with Ncoii scali ng 12811 while initial -distribution is obtained by the solution of classical 
Yang-Mills equations 113ill : also in this case we assume longitudinal boost invariance and y = q. 

In Fig. [T] we plot the time evolution of the particle number for the several initializations we 
consider in this work. On the horizontal axis we measure the time difference r - tq where tq 
















Figure 1: Early time evolution of particle number for thermal Glauber (indigo) and CYM (orange) initializations. 
Initial gluon number corresponds to dNgjdy =1040. 



Figure 2; Time evolution of particle spectra for thermal Glauber (indigo) and CYM (orange) initializations. In both 
simulations the initial gluon number corresponds to dNgjdy =1040. 


corresponds to the initialization time; in this study we ehose tq = 0.6 fm/e for the ease of the 
Glauber initialization and tq = 0.2 fm/e for the CYM initialization. In the latter ease we assume a 
smaller initialization time beeause we do not need to assume any pre-thermalization, whieh instead 
is assumed in the Glauber initialization. In the figure, Nquarks corresponds to the total number of 
quarks plus the total number of antiquarks. The results we present in this article correspond to 
simulations run for Au-Au eollisions at RHIC energy, for b = 2.5 fm; the initial gluon multiplieity 
is dNgjdy =1040. 


5 









Figure 3: Left panel. Early time evolution of gluons and fermions number for the case of CYM initialization and 
several choices for the coupling. Right panel, ratio of total quark number over gluon number as a function of time. 
Data style convention in the right panel is the same as the one in the left panel. Calculations correspond to Au-Au 
collision dX ^fs - 200A GeV and b - 2.5 fm. 


In Fig. [2] we plot the results for initial and final partiele speetra at midrapidity |y| < 0.5. In 
the figure we do not show initial speetra for quarks beeause we assume that at t = Tq there are no 
quarks in the system. This might be a crude approximation for the Glauber initialization since it is 
equivalent to assume that some dynamics took place up to r = Tq to thermafi z e the system in the 
transverse plane, nevertheless affecting only gluons and not producing any quark; nevertheless for 
the early times initializations this might be not a so bad approximation since we are only assuming 
that within a time of the order of 1/2 ^ « Tq the system is mainly made of gluons, then as soon 
as we initialize the dynamics the QCD processes cause the conversion of gluons to quarks. We 
remark that a more realistic condition including a finite number of quarks can only make stronger 
our final result that quark production is quite fast and abundant. 

The results summarized in Figg. [T] and |2] are obtained assuming a value of as = 0.7. As 
a reference we remind that this value of the coupling corresponds to a value of ij/s « 2/An at 
temperature T = 0.3 GeV, which is not far from the value of 77/5 which is normally used in 
ultrarelativistic heavy ion collisions simulations to quantify the dissipation of the QGP fluid. 

In the left panel of Fig. |3]we plot the gluon and total quark numbers as a function of time, 
focusing on the CYM initialization and on the early stage of the time evolution of the fireball. For 
other initializations we obtain similar results which are not shown for space sake. In the figure 
we have shown three different cases corresponding to three different choices for a,. In particular, 
diamonds and circles correspond to a fixed value of as, namely a, = 0.7 and as = 0.3 respectively; 
on the other hand squares correspond to the case of the temperature dependent as given by Eq. ® . 

The definition of a quark formation time Tq is quite arbitrary, therefore we define it as the 
time necessary to have a mixture made of 50% of quarks plus antiquarks and 50% of gluons, 
corresponding to R = I with R = (Nq + Nq)/Ng. As already discussed in the case as = 0.7 the 
system is very efficient in converting gluons to quarks, in fact the total number of quarks becomes 
comparable to the one of gluons within Tq ~ 0.15 fm/c. It is interesting that in the case of the 
perturbative running coupling, indigo data in Fig. |3l even if the average coupling is smaller than 
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Figure 4: Left panel: Early time evolution of gluons and fermions number for the case of CYM initialization and 
several choices for the coupling, for the case of massive quasiparticles. Right panel: ratio of total quark number over 
gluon number as a function of time. Data style convention in the right panel is the same as the one in the left panel. 
Calculations correspond to Au-Au collision si ^fs - 200A GeV and b - 2.5 fm. 


the case as = 0.7 (of about a factor of 1.5 when averaged on the fireball) the system is still quite 
efficient in producing quarks via the QCD inelastic processes. In fact in this case we find Tq « 0.6 
fm/c, which is about a factor of 4 larger than the time required in the case as = 0.7 but still 
very fast compared to the time evolution of the fireball. Finally in the case as = 0.3 we find 
that quark production is quite slow compared to the previous cases, having Tq ~ 2 fm/c which 
is approximately 3.3 times larger than Tq obtained within the perturbative running coupling and 
approximately 1 /2 of the lifetime of the QGP in the fireball. 

In the right panel of Fig. [3] we plot the ratio R = (Nq + Nq)/Ng again for the case of the CYM 
initialization, for the same values of considered in the left panel of the same figure. In the 
case of as = 0.7 the ratio R needs a short time range of about 1 fm/c to reach a value close to 
the equilibrium value for a massles gas, R = 9/4. On the other hand in the case of a running 
perturbative coupling, even if the quark production has been efficient in the early stages after the 
collision, the system is not able to chemically equilibrate onto the chemical equilibrium value 
within the lifetime of the QGP in the fireball; in fact R saturates to a somewhat smaller value, 
7? « 1.7 while 7? = 9/4 at equilibrium. Finally in the case of weak coupling as = 0.3 we find the 
value R ~ 1.2 is reached within the lifetime of the fireball, meaning the system is unable to reach 
chemical equilibration in weak coupling regime even though quark fraction is still abundant. The 
different asymptotic behaviours in this case of massles quasiparticles are due to the freezout which 
freezes the interactions hence arresting conversion processes from quarks to gluons. 


3.2. The massive case: temperature independent masses 

In the left panel of Fig. 0] we plot the gluon and total quark numbers as a function of time, for the 
case of massive quasiparticles. In order to emphasize the role of a finite mass for the quasiparticles, 
as well as of the different mass among quarks and gluons, we have restricted ourselves to the case 
of temperature independent masses. According to the results of 11251] we have chosen the ratio of 
the gluon to quark mass to be 3/2, neglecting a mild temperature dependence of the masses in the 
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Figure 5: Left panel. Early time evolution of gluons and fermions number for the case of CYM initialization and 
several choices for the coupling, for the case of F-dependent masses. Right panel, ratio of total quark number over 
gluon number as a function of time. Data style convention in the right panel is the same as the one in the left panel. 
Calculations correspond to Au-Au collision s.t ^fs - 200A GeV and b - 2.5 fm. 


range of temperature relevant for the evolution of the fireball at RHIC energy. We have ehosen 
the quark mass to be = 0.5 GeV in agreement with the results of OSll . We find that regardless 
of the ehoiee for the eoupling, eonsidering massive quasiparticles amounts to increase the quarks 
to gluons number ratio R of about 30% along with a general acceleration of Tg, compare Figg. [5] 
andlH This is easily understood because in this case quarks are lighter than gluons and statistically 
it is easier to convert a pair of gluons into a qq pair. 

3.3. The massive case: temperature dependent masses 

In the left panel of Fig. [5] we plot the gluon and total quark numbers as a function of time, for 
the case of massive quasiparticles with a temperature dependent mass. The quasiparticle masses 
are computed according to 11251] in order to fit the lattice thermodynamics data; in other words, the 
masses we consider in this section are such that the equation of state of the fluid we simulate is 
that of lattice QCD. Left panel of Fig. [5] focuses on the very early history of the fireball; we find 
that this case does not differ in a relevant way from the case of temperature independent masses 
discussed in the previous subsection. In fact also in this case perturbative QCD processes lead to 
fast quarks production. However a considerable quantitative difference is found for the final values 
of the ratio R of quarks to gluons, which is more evident in the case of strong coupling or^ = 0.7. 
Indeed in the case of temperature independent masses of Fig.|4]we found that R approaches « 2.4 
within few fm/c; instead in the case of temperature dependent masses we find that R « 3 at the end 
of the evolution. 

On the other hand, for the running coupling constant case, circles in Fig. |5l we find that the 
use of the temperature dependent masses does not affect very much the result obtained within 
fixed masses. This can be understood a posteriori since the quark and gluon masses in Fig. |4] 
do not differ too much from those of Fig. [51 besides the peripheral cells where the temperature 
is smaller thus the masses are larger, but the coupling is smaller than the case as = 0.7 (roughly 
corresponding to q!s = 0.15) implying the system is less effective in converting gluons to quarks. 
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Figure 6: Comparison of the ratio R for the Glauber and the CYM initializations. Left panel corresponds to Glauber 
initialization, right panel to CYM initialzation. Calculations are performed with a fixed as - 0.7, for Au-Au collision 
at = 200A GeV and b - 2.5 fm. 



Figure 7: Comparison of the ratio R for the Glauber and the CYM initializations. Left panel corresponds to Glauber 
initialization, right panel to CYM initialzation. Calculations are performed with a temperature dependent as, for 
Au-Au collision at V* = 200A GeV and b - 2.5 fm. 


The regions of the fireball where temperature is elose to the eritieal temperature are still important, 
beeause the smaller the temperature the larger the ratio R in this ease, but these regions ean be 
elFeetive in the produetion of quarks only if the eoupling is large enough, whieh is aehieved by the 
ease = 0.7 but non in the ease of the running perturbative eoupling. This explains the results on 
the right panel of Fig. [51 

3.4. CYM versus Glauber initializations 

In this subseetion we eompare the evolution of the quarks to gluons ratio obtained by the 
Glauber and the CYM initializations, for both oases of strong fixed and temperature dependent 
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couplings. The cases of both temperature dependent and constant masses are considered. This 
comparison is instructive since it permits to enlighten about the role of the different initialization 
on quark production. It is useful to remind here that the case of temperature dependent masses 
corresponds to assume the QGP fluid is characterized by the equation of state of lattice QCD. 

In Fig. 0 we plot the evolution of 7? as a function of time for the Glauber initialization (left 
panel) and the GYM initialization (right panel). The results are shown for the cases of constant 
mass and temperature dependent mass, with perturbative running and fixed strong coupling (the 
latter corresponding to as = 0.7). In the case of strong coupling, left panel of Fig. we find that 
for the Glauber initialization switching from a constant mass to a thermal one leads to an increase 
of the ratio 7? at later times. This is mainly due to the fact that the number of gluons decreases 
(and that of quarks increases) continuously in the case of a temperature dependent mass, while it 
saturates after a time of about 2 fm/c in the case of a constant mass. We find the same result for 
the GYM initialization, even if in this case the effect of turning on a temperature dependent mass 
is milder. However, the conclusion is that regardless the initial condition chosen in the simulation, 
at later times the fireball consists of about 80% of quarks and only 20% of gluons even starting 
from a condition in which the system is made of gluons only. Finally in Fig.|7]we plot the ratio 
for the case of as(T). 

4. Conclusions 

In this Article we have studied quark production from a pure gluon plasma produced in ul- 
trarelativistic heavy ion collisions. Our main goal has been to study how a pure gluon plasma 
evolves dynamically to a quark-gluon plasma by means of QCD inelastic processes; the impact 
of the initial condition on the relative composition in quarks and gluons and on quark production 
times has been investigated. We have considered two different kinds of initializations: the usual 
Glauber initialization and one based on the gluon spectrum obtained by the classical Yang-Mills 
theory (CYM) of the glasma. Given the initial condition, we have evolved it according to relativis¬ 
tic transport theory, and quark production has been taken into account by means of QCD inelastic 
processes. 

The definition of a quark formation time is quite arbitrary, therefore we have defined it as 
the time necessary to have a mixture made of 50% of quarks plus antiquarks and 50% of gluons. 
We have found that the initial condition does not have a relevant impact on this is true both 
if we consider a pure perturbative coupling and a strong coupling regime in which as is chosen 
scale-independent and of magnitude such that the corresponding rj/s is in the range 1-^2. The 
coupling has some effect on r^: changing the coupling from strong to weak we have obtained 
Tg increases about of a factor 4; nevertheless even in weak coupling we find abundant fermion 
production takes place within r < 1 fm/c. 

We have then computed the ratio = (Ng+Ng)/Ng and we have found it is affected not only by 
the coupling strength but also on the kind of quasiparticle mass we implement in the simulation. 
In the case of strong coupling we have found that regardless of the initial condition the ratio R 
approaches about 2.5 3.5; in the case of weak coupling the effect of the initialization is quite 

strong for the case of temperature dependent masses, being instead not so strong in the case of 
massive quasiparticles with a constant mass. It is interesting that even in the weak coupling case 


10 


the asymptotic R x 2 ^ 2.5 signaling an abundant production of quarks plus antiquarks has been 
achieved during the fireball’s lifetime. 

Among the many interesting aspects which deserve further study we briefly mention here the 
possibility to include, in future studies, the Bose-Einstein enhancing factors in the collision integral 


which might lead to a transient condensate in the very early stages of the glasma evolution [|32L I33 


3J,l35Ll3^, hence affect the conversion of gluon to quarks 113711 . Moreover it will be interesting to 


compare the quark production in QGP studied in the present article with the one obtained by the 
shattering of the flux tubes via the Schwinger mechanism, following the lines of 113811 . Both these 
problems will be the subject of forthcoming publications. 
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